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We report a search for charmless hadronic decays of neutral B mesons to the final state K'^ K^-k*^ . 
The results are based on a 711 fb^^ data sample that contains 772 x 10® BB pairs, and was collected 
at the T(4S) resonance with the Belle detector at the KEKB asymmetric-energy e"'"e~ collider. We 
find the first evidence for this decay with a significance of 3.5 standard deviations and measure its 
branching fraction as B{B° K+K'-k^) = [2.17 ± 0.60(stat) ± 0.24(syst)] x 10"^ 

PACS numbers: 13.25.Hw, 14.40. Nd 



The i3-meson decay — >■ K^K^tt^ is suppressed in 
the standard model (SM) and thus offers a useful probe 
for new physics beyond the SM. Figure [1] shows typical 
Fcynman diagrams that contribute to this decay. The 
dominant one is the color- and Cabibbo-suppressed h ^ u 
tree transition followed by the internal W exchange dia- 
gram leading to 5° K*^K^ with K*^ K^Tr°. The 
latter diagram dominates in the decay — > K^K^ , for 
which only upper limits have been placed on the branch- 
ing fraction \^]^ . This is in contrast to the related decays 
(having two kaons in the final state) that are already ob- 
served such as B° K°K°, B+ K°K+ 0,11], and 



where the b d gluonic penguin 
amplitude can contribute as well Q . 

The three-body decay B^ — > K^K'-k^ has not yet 
been observed, with only one measured upper limit of 
B{B^ /v+A'-tt") < 19 X 10-6 at 90% confidence level 
from the CLEO Collaboration Q . Intermediate resonant 
modes that decay preferentially to this final state have 
also not been seen. A search for a related channel by Belle 
has set an upper limit of B{B° (l>Tr°) < 1.5 x 10"^ 0. 
The latter mode is quite sensitive to possible beyond- 
the-SM contributions; a branching fraction of ©(lO"'') 
would constitute evidence for new physics |ll| . No exper- 
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FIG. 1: Typical Feynman diagrams that contribute to the 
decay K^K~n'^: (a) 6 — >■ w tree and (b) internal W 

exchange. 



imental information is available for other potential reso- 
nance modes such as K*{892)^K'^, K^{U30)^K'^, and 
/o(980)7rO. Pqj. ^-i^e ^^^^^ qo _^ K*{m2)^K^ dominated 
by internal W exchange [Fig. [Ijb)], the branching frac- 
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tion is predicted to be in the range 10~^ to 10~^ 

Another motivation for the study of — ;> K^K^tt 
comes from the observation of — ^ K^K^n^ by 
the BaBAR Collaboration In particular, an unex- 
pected structure is seen near 1.5GeV/c^ in the K'^K~^ 
invariant-mass spectrum, which accounts for about half 
of the total events. Similar structures have also been 
observed in the Dalitz plots of i?+ — > K^K~K^ and 
B^ — !> K^K~K^ decays [l5-17|. If these structures are 
due to a particular K~ resonant state, it should show 
up in B^ — > K^K^TT^] on the other hand, if it is a re- 
flection from the b d penguin, it will not contribute 
to A'+A'^tt". Since the u and d quarks are spectators in 
the b ^ u tree diagram [Fig. [Ija)] for B^ — !> A'+A'^7r+ 
and —> A'+A'~7r°, respectively, one can estimate the 
branching fraction for the latter using the BaBAR results. 
Assuming isospin symmetry and the b ^ u transition to 
be the main contributor to B^ — >■ AT+AT^tt^, we expect 
its branching fraction to be at the level of 3 x 10~^, which 
is well within Belle's reach. 

Our results are based on a data sample containing 
772 X 10^ BB pairs collected at the Ti^AS) resonance with 
the Belle detector [3| at the KEKB asymmetric-energy 
e+e" (3.5 on 8.0 GeV) collider The principal detec- 
tor components used in the study are a silicon vertex de- 
tector, a 50- layer central drift chamber (CDC), an array 
of aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and a CsI(Tl) crystal electromagnetic calorimeter 
(ECL). All these components are located inside a 1.5 T 
solcnoidal magnetic field. 

To reconstruct B^ — !> A'+A'^tt" decay candidates, we 
combine two oppositely charged kaons with a tt" meson. 
Each track candidate must have a minimum transverse 
momentum of lOOMeV/c, and a distance of closest ap- 
proach with respect to the interaction point of less than 
0.2 cm in the transverse r-(j) plane and less than 5.0 cm 
along the z axis, where the z axis is defined by the di- 
rection opposite the e+ beam. Identification of charged 
kaons is based on a likelihood ratio Rk/t^ ~ Ck+c ' 



where Ck and C-n denote the individual likelihoods for 
kaons and pions, respectively, calculated using specific 
ionization in the CDC, time-of-fiight information from 
the TOF, and the number of photoelectrons from the 
ACC. A requirement Rk/-k > 0.6 is applied to select 
both kaon candidates. The kaon identification efficiency 
is approximately 86% and the probability of misidenti- 
fying a pion as a kaon is 11%. We reconstruct 7r° can- 
didates from photon pairs that have an invariant mass 
between 112 and 156MeV/c^, corresponding to ±3.5cr 
around the nominal tt*^ mass (20j . These photons are 
reconstructed from neutral clusters in the ECL with en- 
ergy above 60 (100) MeV in the barrel (endcap) region. 
In addition, requirements on the tt'^ decay helicity angle, 
|cos6'hGi| < 0.95, and the 7r° mass-constrained fit statis- 
tic, Xmass < 50, are imposed. Here, 6'hei is the angle 
between one of the daughter photons and the B momen- 
tum in the TT*^ rest frame. 

i?-meson candidates are identified using two kine- 
matic variables: beam-energy constrained mass, Mhc ~ 

and energy difference, AE — 



Ey, 



where -Bbeam is the beam energy, and 



Pi and Ei are the momentum and energy, respectively, of 
the ith daughter of the reconstructed B in the center-of- 
mass (CM) frame. We retain events with 5.271 GeV/c^ < 
Afbc < 5.289 GeV/c2 and -0.30 GcV < A£; < 0.15 GeV 
for further analysis. The Mhc requirement corresponds 
to approximately ±3cr around the nominal B'^ mass [20j : 
we apply a looser window of (—12(7, +6a) around AE = 
because it is used in the fitter (as described below). The 
average number of B candidates found per event is 1.3. In 
events with multiple B candidates, we choose the one(s) 
whose 7r° has the lowest Xmass value. If more than one 
B candidate shares the same 7r° meson, the candidate 
yielding the best vertex fit is selected. 



The dominant background is from the e+e~ qq 
{q = u, d, s, c) continuum process. To suppress this back- 
ground, observables based on the event topology are uti- 
Hzed. The event shape in the CM frame is more spherical 
for BB events and jetlike for continuum events. We em- 
ploy a neural network [2lj to combine the following six 
input variables: the Fisher discriminant formed from 16 
modified Fox- Wolfram moments 221, the cosine of the 



angle between the B momentum and the z axis, the co- 
sine of the angle between the B thrust and the z axis, the 
cosine of the angle between the thrust axis of the B can- 
didate and that of the rest of the event, the ratio of the 
second- to the zeroth-order Fox- Wolfram moments (all 
of these quantities being calculated in the CM frame), 
and the separation along the z axis between the vertex 
of the B candidate and that of the remaining tracks. The 
training and optimization of the neural network are ac- 
complished with signal and qq Monte Carlo (MC) simu- 
lated events. The signal MC sample is generated with the 
EvtCen program [23[ by assuming a three-body phase 
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space. We require the neural network output (Cnb) to 
be above 0.2 to substantially reduce the continuum back- 
ground. The relative signal efficiency due to this require- 
ment is approximately 88%, whereas the continuum sup- 
pression achieved is close to 92%. The remainder of the 
Cnb distribution peaks strongly near 1.0 for signal, and 
thus we have difhculty in modeling it with an analytic 
function. However, its transformed variable 



an irreducible correlation with AE owing to shower leak- 
age in the ECL. The extended likelihood function is 



£ = exp - ^ X ]J 



(3) 
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where CMB,-min = 0.2 and Cnb,) 
tion with a Gaussian-like tail. 

The background due to B decays via the dominant 
b ^ c transition is studied with an MC sample of a 
collection of such decays. The resulting Mbc distribu- 
tion is found to peak strongly in the signal region. We 
also observe two peaks in the K~^K~ invariant-mass 
spectrum that corresponds to the contributions from (a) 
D'^ K^K~ peaking at the nominal mass [i^], and 
(b) — > K~7T'^ with the peak shifted slightly from the 
D'^ mass owing to K-n misidentification. To suppress 
these peaking contributions, we exclude candidates for 
which the invariant mass of the K~^K~ system lies in 
the range of [1846, 1884] MeV/c^ (about ±5cr around the 
nominal 1?° mass). In the case of (b), we use the pion 
hypothesis for one of the tracks. The surviving events 
constitute the "generic BB" background. 

There are a few background modes that contribute 
in the Mbc signal region having the AE peak shifted 
to positive values. The so-called "rare peaking" back- 
ground modes arising mostly from K~tt misidentification 
are identified with a BB MC sample in which one of 
the B mesons decays via b — > u,d, s transitions with 
known or estimated branching fractions. The rare peak- 
ing background includes the K~^Tr~Tr'^ nonres- 
onant decay as well as possible intermediate resonant 
modes that result in the if+7r~7r° final state, such as 
B" A'*(892)07rO and B° ii:*(892)+7r-. The events 
that remain after removing the signal and rare peak- 
ing components comprise the "rare combinatorial" back- 
ground. 

The signal yield is obtained with an unbinned extended 
maximum likelihood fit to the two-dimensional distribu- 
tions of AE and C'j^g. We define a probability density 
function (PDF) for each event category j (signal, qq, 
generic BB, rare peaking, and rare combinatorial BB 
backgrounds): 



where rij is the yield of event category j. The correctly 
reconstructed (CR) and misreconstructed fragments of 
the i3-meson decay referred to as self-crossfeed (SCF) 
components of the signal are considered distinct in the 
fitter: their combined PDF is n^igX [(1-/) Pcr+/Pscf], 
= 1.0, has a distribu- where ngjg is the total signal yield and / is the SCF frac- 



Cnb 



(1) 



P]^r,iAE^)V,{C^s), 



(2) 



where i denotes the event index. Since the correlation be- 
tween AE and C'j^g is found to be negligible, the product 
of two individual PDFs is a good approximation for the 
combined PDF. We apply a tight requirement on Alhc 
rather than including it in the fitter because it exhibits 



tion fixed to the MC expected value of 3%. 



TABLE I: List of PDFs used to model the AE and C'nb dis- 
tributions for various event categories. G, AG, CB, and PoIy2 
denote Gaussian, asymmetric Gaussian, Crystal Ball [24| . and 
second-order Chebyshev polynomial function, respectively. 



Event category 



AE 



CR signal 
SCF signal 
Continuum qq 
Generic BB 
Rare peaking BB 



CB+AG 
histogram 
Poly2 
Poly2 
2G 



SAG 
histogram 
AG 
AG 
AG 



Rare combinatorial BB histogram 3 AG 

Table H] lists the PDF shapes used to model the AE 
and C'j^^ distributions for each event category. Distri- 
butions that are difficult to parametrize analytically are 
modeled with histograms. The yields for all event cat- 
egories except the rare peaking BB background are al- 
lowed to vary in the fit. We fix the yield of the rare 
peaking BB component to the value calculated using 
the branching fraction measured in an amplitude anal- 
ysis of B° K+Tr-n° The following PDF shape 



parameters of the qq background arc floated: the two pa- 
rameters of the second-order Chebyshev polynomial used 
for AE, and the mean and two widths of the asymmet- 
ric Gaussian function used to model C^^. The PDF 
shapes for signal and other background components are 
fixed to the corresponding MC expectations. We adjust 
the parameters of the signal AE and C'pfg PDFs to ac- 
count for possible data-MC differences, according to the 
values obtained with a large-statistics control sample of 
_B+ D"(A'+7r~7r°)7r"'". The same correction factors are 
also applied for the rare peaking BB background. 

Figure [2] shows the AE and C'^g projections of the fit 
applied to 39066 candidate events. We obtain 299 ± 83 
signal events («sig), 32167±428 continuum qq, 3814±517 
generic BB, and 2691 ±321 rare combinatorial BB 
background events. The statistical significance of the 
signal is 3.8 standard deviations. It is calculated as 
■\/2 log(£o/£max), where £o and £max are the fit likeli- 
hood values with the signal yield set to zero and the best- 
fit case, respectively. The obtained background yields are 
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FIG. 2: (color online). Projections of candidate events onto 
(left) AE for C'mb > 3 and (right) C^vs for \AE\ < 30MeV. 
Points with error bars are the data, solid (blue) curves are the 
total PDF, dashed (red) curves are the total background, dot- 
ted (green) curves are the sum of continuum qq and generic 
BB backgrounds, dash-dotted (magenta) curves are the con- 
tinuum qq background, and filled (cyan) regions show the sig- 
nal. 



consistent with the respective MC predictions. The sig- 
nal decay branching fraction is calculated as 



(4) 



where N^g is the total number of BB pairs (772 x 10^), 
Ercc is the signal reconstruction efficiency (19.6%) ob- 



tained in the study described below, and 



K/ir 



denotes 



the kaon-identification efficiency correction factor that 
accounts for a small data-MC difference. It is given by 



— -data 



MC 



(5) 



where ej^j!^ (^A/ir) ^^'^ efficiency of the Rk/tt require- 
ment in data (MC simulations). The rj^/^ value per kaon 
track is 0.95, resulting in a total t^/tt = 0.95^ = 0.90 for 



two kaons. We have verified that the R 



correction 



factor is almost constant over the Dalitz plot. For the 
branching fraction calculation presented in Eq. (j4|), we 
assume equal production of B°B^ and B^B^ pairs at 
the T(4iS') resonance. The resulting value is 



B{B" 



K+K 



[2.17 ± 0.60 ± 0.24] X 10-^ (6) 



where the uncertainties are statistical and systematic, re- 
spectively. The contributions to the systematic uncer- 
tainty are discussed below and listed in Table III 

The uncertainties due to the PDF shape parameters 
are estimated by varying all fixed parameters by ilcr. 
To assign a systematic error for the histogram PDF used 
to model Ai? for the rare combinatorial component, we 
carry out a series of fits by fluctuating each of the his- 
togram bin contents according to the Poisson distribu- 
tion. The spread of the fitted signal yields is taken as the 
systematic error. We also vary the yield of final states 
that dominantly contribute to that component according 



to their errors. As we use a fairly complex function (a 
sum of three asymmetric Gaussians) to model the sig- 
nal C'ffg PDF shape, we evaluate possible systematics 
due to the uncertainty in the functional dependence by 
checking other alternatives. This systematic contribution 
is denoted as "Signal C'j^^ functional dependence" in Ta- 
blc|TTl The uncertainty due to the fixed (small) SCF frac- 
tion is estimated without knowing a priori how these SCF 
events vary across the Dalitz plot. We adopt a conserva- 
tive approach to vary the SCF fraction by ±50% when 
calculating the associated systematic error. The poten- 
tial fit bias is evaluated by performing an ensemble test 
comprising 200 pseudocxperiments, where the signal and 
rare peaking background components are embedded from 
the corresponding MC samples, and the PDF shapes are 
used to generate the data for the other event categories. 
We obtain an almost Gaussian pull distribution of unit 
width, and add the mean and error on the pull in quadra- 
ture for assigning the systematics. Uncertainty due to 
continuum suppression is derived with the control sample 
by comparing the nominal fit result with that obtained 
without any Cnb requirement. We estimate the error due 
to the Mbc requirement by varying its nominal selection 
threshold by the resolution. The D*+ D°{K- t:+)tt+ 
control sample is used to determine the systematic un- 
certainty due to the Rk/tt requirement. The systematic 
uncertainty due to 7r° reconstruction is evaluated by com- 
paring data-MC differences of the yield ratio between 
77 — 7r°7r°7r" and 77 — 7r+7r~7r°. We use partially recon- 
structed D*^ £'"(/\g7r+7r^)7r+ decays to assign the 
systematic uncertainty due to cliarged-track reconstruc- 
tion (0.35% per track). To account for the possible vari- 
ation of efficiency across the Dalitz-plot distribution, we 
calculate a weighted signal reconstruction efficiency by 
fitting different regions of that distribution. The mean 
value is used to obtain the branching fraction and the 
error is taken as the systematic contribution due to the 
efficiency variation. The total systematic uncertainty is 
calculated by summing all these uncertainties in quadra- 
ture. To determine the significance of our measurement, 
we use a convolution of the statistical likelihood with 
a Gaussian function of width equal to the additive sys- 
tematic errors that only affect the signal yield. The total 
significance, including these uncertainties, is 3.5 standard 
deviations. 

To elucidate the nature of the observed signal, es- 
pecially whether there are contributions from the de- 
cays with intermediate resonant states, we study the 
K'^K~ and K'^t:^ invariant mass distributions. We per- 
form the [AS, C'j^g\ two-dimensional fit in bins of the 
m{K^ K^) and in{K^Tr^) distributions after applying 
the orthogonal requirements m{K^'K'^) > 1.5GeV/c^ and 
■m{K'^K^) > 2.0GeV/c^, respectively. These require- 
ments suppress kinematic refiections. Figure [3] shows the 
resulting signal yields along with their statistical errors. 
With these data, we cannot make any definitive state- 
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TABLE II: Summary of various systematic uncertainties. Tfie 
first and second fiorizontal blocks denote tlie additive and 
multiplicative systematic uncertainties, respectively. 



Source 


Uncertainties (%) 


Signal PDF 


+3.4 


-2.9 


Generic B B PDF 


+2.4 


-3.1 


Combinatorial background PDF 


+ 1.3 


-2.0 


Peaking background PDFs 


+ 1.7 


— 1.9 


Fixed liistogram PDF 


+1.7 


-2.0 


Signal Cjvjj functional dependence 


+2.3 


-2.3 


Fixed SCF fraction 


+ 1 7 


— 1.7 


Fit bias 


+2.4 


-2.4 


Continuum suppression 


+2.2 


-2.2 


Requirement on Afbc 


+ 1.5 


-0.2 


Kaon ID requirement 


+ 1.9 


-1.9 


tt" detection efficiency 


+4.0 


-4.0 


Charged track reconstruction 


+0.7 


-0.7 


Efficiency variation over Dalitz plot 


+7.5 


-7.5 


Number of BB pairs 


+ 1.4 


-1.4 


Total 


+ 11.1 


-11.3 



ment about possible intermediate K~ resonances, in- 
cluding the structure seen by BaBAR near 1.5 GeV/c^ Q. 
It is worth noting here that the recent LHCb study of 
B± -> K+K-n^ decays has revealed an unidenti- 
fied structure in the same mass range; however, it is only 
present in events, giving rise to a large local CP asym- 
metry. Furthermore, we observe some excess of events 
around 1.4GeV/c^ in the A'+tt" invariant-mass spectrum. 
A detailed interpretation will require an amplitude anal- 
ysis with higher statistics that would be available at a 
next-generation flavor factory [2^ . 
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dence for the decay. 
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FIG. 3: Signal yield distributions as a function of (left) 
m{K+K-) with m{K+TT°) > l.SGeV/c^ and (right) 
miK+Tv") with m{K+K-) > 2.0GeV/c2. Each point is ob- 
tained from a two-dimensional [AE,C'pjg] fit. 

In summary, we report measurement of the suppressed 
decay ^ K+K'tt^ using the full T{AS) data sam- 
ple collected with the Belle detector. We employ a two- 
dimensional fit for extracting the signal yield. Our mea- 
sured branching fraction B{B^ K^K'tt^) ~ [2.17 ± 
0.60(stat) ± 0.24(syst)] x 10"^ constitutes the first evi- 
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